home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dgttrf.z / dgttrf
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))                                                          DDDDGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DGTTRF - compute an LU factorization of a real tridiagonal matrix A using
  10.      elimination with partial pivoting and row interchanges
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DGTTRF( N, DL, D, DU, DU2, IPIV, INFO )
  14.  
  15.          INTEGER        INFO, N
  16.  
  17.          INTEGER        IPIV( * )
  18.  
  19.          DOUBLE         PRECISION D( * ), DL( * ), DU( * ), DU2( * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      DGTTRF computes an LU factorization of a real tridiagonal matrix A using
  23.      elimination with partial pivoting and row interchanges.
  24.  
  25.      The factorization has the form
  26.         A = L * U
  27.      where L is a product of permutation and unit lower bidiagonal matrices
  28.      and U is upper triangular with nonzeros in only the main diagonal and
  29.      first two superdiagonals.
  30.  
  31.  
  32. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      DL      (input/output) DOUBLE PRECISION array, dimension (N-1)
  37.              On entry, DL must contain the (n-1) subdiagonal elements of A.
  38.              On exit, DL is overwritten by the (n-1) multipliers that define
  39.              the matrix L from the LU factorization of A.
  40.  
  41.      D       (input/output) DOUBLE PRECISION array, dimension (N)
  42.              On entry, D must contain the diagonal elements of A.  On exit, D
  43.              is overwritten by the n diagonal elements of the upper triangular
  44.              matrix U from the LU factorization of A.
  45.  
  46.      DU      (input/output) DOUBLE PRECISION array, dimension (N-1)
  47.              On entry, DU must contain the (n-1) superdiagonal elements of A.
  48.              On exit, DU is overwritten by the (n-1) elements of the first
  49.              superdiagonal of U.
  50.  
  51.      DU2     (output) DOUBLE PRECISION array, dimension (N-2)
  52.              On exit, DU2 is overwritten by the (n-2) elements of the second
  53.              superdiagonal of U.
  54.  
  55.      IPIV    (output) INTEGER array, dimension (N)
  56.              The pivot indices; for 1 <= i <= n, row i of the matrix was
  57.              interchanged with row IPIV(i).  IPIV(i) will always be either i
  58.              or i+1; IPIV(i) = i indicates a row interchange was not required.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))                                                          DDDDGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      INFO    (output) INTEGER
  75.              = 0:  successful exit
  76.              < 0:  if INFO = -i, the i-th argument had an illegal value
  77.              > 0:  if INFO = i, U(i,i) is exactly zero. The factorization has
  78.              been completed, but the factor U is exactly singular, and
  79.              division by zero will occur if it is used to solve a system of
  80.              equations.
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.